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ABSTRACT 

The  model  for  occurrence  of  fingerprint  characteristics  in  terms 
of  multinomial  trials  on  a grid  of  cells  is  extended  to  consider 
dependence  between  the  cells.  The  occurrence  of  the  characteristics 
is  modelled  as  a two-dimensional  Markov  process. 
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AS  A TWO-DIMENSIONAL  PROCESS 

Stanley  L.  Sclove 
University  »f  Illinois  at  Chicago  Circle 


1.  INTRODUCTION 


The  individuality  of  a fingerprint  is  based  on  the  pattern  of  occurrence 
of  the  ridge-line  details.  These  minutiae,  called  Galton  [2]  characteristics,  are 
of  ten  types:  islands,  bridges,  spurs,  dots,  ridge  endings,  forks  (bifurcations), 
lakes  trifurcations,  double  bifurcations,  and  deltas.  (See  [3]  for  diagrams  and 
detailed  descriptions.)  In  [3]  the  estimation  of  fingerprint  probabilities  based 
on  Galton  characteristics  was  treated  according  to  the  following  model: 

Assimption  1.  A fingerprint  is  considered  in  terms  of  a grid  of  one  millimeter  cells. 
Assumption  2.  For  each  cell  of  the  grid  there  are  13  possibilities:  either  the 
cell  is  empty,  or  one  of  the  following  12  possibilities  has  occurred:  island,  bridge, 
spur,  dot,  ending  ridge,  fork,  lake,  trifiircation,  double  bifurcation,  broken  ridge 
(two  ridge  endings),  or  some  other  multiple  occurrence. 

Assumption  3.  There  is  statistical  independence  between  cells. 

Under  this  model  the  probability  distribution  for  a given  cell  is  the  point 
multinanial 


^0  ’^l  ^12 

^0  ^1  ‘“^12 

where,  for 

i = 0 , 1 } 2 , • • • ) 12, 

or  not , so 

Zq  + Zi  + ...  + z^2 

the  probability  P of  a given  configuration  of  t cells  is  the  product  over  cells, 

p - n ^12, c / V 

p = n Pq  p^  ...p^2  » (1-1) 


where,  for  i = 0,1,2, ... ,12,  = 1 or  0 according  as  the  i-th  possibility  occurs 

in  the  c-th  cell  or  not. 
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Thus  the  probability  for  a conf igu-ation  with  empty  cells,  cells 
containing  an  island,  k^  cells  containing  a bridge,  ...,  k^^  cells  containing  a delta, 
k^^  cells  containing  two  ridge  endings,  and  k^^  cells  containing  some  other  multiple 
occurrence  is 

ko  k^  k^2 

P = Pq  Pi  ...p^2  ’ 

where,  for  i = 0,1,2, ... ,12,  k^  = number  of  occurrences  of  possibility  i 

across  the  t cells.  The  estimates  of  the  p^'s  (from  [3])  are  given  in  Table  1. 

[INSERT  TABLE  1 HERE.] 

For  example,  consider  the  configuration  of  Figure  1.  It  has  1*3  cells,  37 
of  them  empty,  the  other  six  being  occupied  by  1*  ridge  endings  and  2 forks.  The 


estimated  probability  is 

P = Pq'^P^P^  = .766^'''. 0832**. 0382^;  -IoE^qP  = 11. »*  . 

[INSERT  FIGURE  1 HERE.] 

The  purpose  of  the  present  paper  is  to  study  the  extent  of  departure  from 
Assumption  3 and  to  refine  the  model  according  to  that  departure.  Accordingly, 
the  occurrence  of  fingerprint  characteristics  is  modelled  as  a two-dimensional  process 
to  take  into  account  the  dependence  between  cells. 

Appendix  A gives  a data-analysis  relating  to  dependence  between  cells. 

It  shows  that  the  probability  that  a cell  is  occupied  increases  monotonically 


with  the  number  of  neighbors  occupied. 


2.  MODELLING  DEPENDENCE  AMONG  CELLS 


Let  the  cells  be  numbered  in  some  fixed  order,  say,  as  one  reads 

English,  starting  with  the  top  row  and  moving  from  left  to  right  within  each  row. 

Let  X be  a random  vector  giving  the  outcome  in  the  c-th  cell, 
c 

^c  * ^^0c’^lc’^2c’'“’^12,c^’  ° 


Then 


P = P(X^=x^)P(X2=X2|x^=x^)P(X^=x^|X2=X2,X^=x^) 

. . . P( X^=x^  I X^^^=x^  • • • jX^— x^ ) . 

Under  Assumption  3,  independence  among  cells,  expression  (2.1)  simplifies  to 
P = P(X^=x^)P(X2=X2)...P(X^=x^), 
which  is  the  same  as  (l.l)  since 


(2.1) 


Zq  z z 

P(^~^)  ~ PI^q > ’ ■ ■ ' ’ ^^2 ^ ~ ^ ^0 ’ ^1’ ' " ” * ^12 ^ ^ ~ ^0  ^1  ^12  * 

As  a step  toward  modelling  dependence,  we  introduce 

Assumption  3'.  The  outcome  in  the  c-th  cell  depends  upon  the  outcomes  in  the  other 
cells  only  through  the  outcomes  in  the  adjacent  cells. 

Due  to  the  fact  that  the  probability  (2.1)  forces  one  to  use  a linear  ordering 

of  the  cells,  one  must  write  things  in  terms  of  the  four  prec eding  adjacent  cells 


rather  than  all  eight  adjacent  cells.  More  precisely,  under  Assumption  3',  the 
conditional  probability  P(X  =x  |x  ,X  ,...,X  ) will  not  depend  upon  all  of 

0 C C ^ 

X^  X^  g,  ...,  X^  but  only  upon  four  of  these  variables,  namely,  those  correspon- 
ding tb  the  cell  to  the  left  (west)  of  cell  c,  the  cell  above  (north  of)  cell  c, 
the  cell  just  northwest  of  cell  c,  and  the  cell  just  northeast  of  cell  c.  [If  the 
configuration  were  rectangular  and  indexed  as  (i,j),  then  the  four  cells  upon  which 
the  outcome  in  cell  (i,j)  would  depend  would  be  cells  (i,j-l),  (i-l,j-l),  (i-l,j), 

and  (i-l,j+l).]  If  W denotes  the  matrix  whose  columns  are  these  four  neighbors 

c 


of  X^,  then  Assumption  3'  is 


(2.2) 


Assumption  3'  may  be  viewed  as  an  assiunption  that  the  process  is  a Markov  process; 
see  Appendix  B. 

2.1.  Data  Analysis 

Sets  of  5 cells  were  examined  to  study  the  dependence  of  X^  on  its  four 
preceding  neighbors.  According  to  the  model,  these  sets  were  of  the  form  ^ ^ ^ 


where  y denotes  the  dependent  cell  and  the  x's  its  four  preceding  neighbors. 
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For  any  such  set  s,  let  the  variable  = 1 or  0 according  as  the  fifth  ("y") 

cell  in  the  s-th  set  is  occupied  or  not,  and  let  be  the  number  of  preceding 

adjacent  cells  ("x"  cells)  which  are  occupied;  is  between  0 and  h.  Table  2 

gives,  for  each  value  of  A,  the  proportion  of'  B's  that  are  equal  to  1.  We  say 

[INSERT  TABLE  2 HERE.] 

that  a cell  has  one  adjacency  for  each  of  the  four  preceding  adjoining  cells 

which  is  occupied.  We  shall  refer  to  the  variable  A as  the  number  of  adjacencies. 

The  probability  of  occupancy  increases  monotonically  with  the  number  of  adjacencies. 

Such  absolute  consistency  was  not  expected,  firstly  because  perfect  consistency  seems 

so  rare  in  data  analyses  and  secondly  because  it  was  thought  that  occurrences  in  most 

of  the  four  adjacent  cells  might  crowd  out  occurrence  in  the  fifth  cell. 

Combining  (2.2)  with  the  multinomial  gives 

12 

P(X  =x  |W  ) = n [p.(W  )]^i  , (2.3) 

c c'  C 1C 

where  z = ( z. . z. , . . . , z, ^) . The  model  we  shall  use  for  p.(W  ) is  that  it  depends 
0 1 12  1 c 

on  W^  only  through  the  number  of  cells  occupied.  That  is, 

P. (W  ) = P. (a), 

1 c ^ 

where 

a = a(W  ) 
c 

is  the  number  of  adjacencies  for  cfell  c — the  number  of  occupied  cells  among 

the  four  cells  preceding  and  adjacent  to  cell  c;  the  quantity  a is  either 

0,1, 2, 3 or  U.  We  now  have 

t 12  z. 

P = 11  n [p.(a  )]  (2.U) 

c=l  i=0  ^ ° 

Some  adjustments  are  necessary  for  border  cells  — cells  in  the  first  row,  first  col\imn, 
or  last  column.  See  the  next  section. 

Let  E be  the  event  that  a given  cell  is  occupied.  We  shall  assume  that 
the  probability  of  Possibility  i in  any  given  cell,  which  was  p^  in  the  model  of  [3], 


r 


is 


Assumption  U:  p^(a)  = p^P(ElA=a)/P(E),  i = 1,2,. ..,12. 


p (a)  = 1 - E p (a) 

^ i=l  ^ 

12 

= 1 - [P(E|A=a)/P(E)]  I p.  . 

i=l  ^ 

Note  that,  for  i = 1,2,..., 12,  p^/P(E)  is  simply  the  conditional  probability  of 

occurrence  of  possibility  i,  given  that  the  cell  is  occupied,  so  that  the  effect 
of  Assumption  h is  simply  to  allocate  P(E|A=a)  to  the  twelve  different  possibilities 
in  the  same  way,  regardless  of  the  value  of  A. 


This  gives 

t 12  z. 

p = n n [p.p(E|A  =a  )/:(E)] 
c=li=0  ^ ° ^ 

U 12  k.  -(k  +k_+. . .+k  ) U m 

= [ n p-(a)  ° ] [ n p,  P(E)  ^ ^ [ n P(EiA=a)  , 


where,  for  a = 0,1,2,3,*<, 

kpCa)  = number  of  empty  cells  with  exactly  a adjacencies, 
for  i = 1,2, ... ,12 , 

k^  = number  of  cells  containing  possibility  i 
(as  above),  and  for  a = 0,1, 2, 3,^, 

m^  = number  of  occupied  cells  with  exactly  a adjacencies. 

2.2.  Border  Cells 

Cells  at  the  border , not  being  touched  by  the  full  complement  of  four 
preceding  adjacent  cells,  require  some  special  treatment.  One  could  take  the  results 
in  border  cells  as  given  and  take  the  probabilities  for  the  other  cells  conditionally 
on  the  outcomes  in  the  border  cells,  but  this  would  result  in  considerable  reduction 
in  the  effective  sample  size.  (E.g.,  l8  of  the  *43  cells  in  Figure  1 are  border  cells.) 


We  wish  to  use  such  information  as  is  pi'csent.  I.e.,  we  shall  make  use  of  the 
fact  that  some  hut  not  all  U preceding  adjacent  cells  are  represented.  E.g., 
if  there  are  3 adjacent  cells  and  2 are  occupied,  then  we  know  that  if  the  U cells 
were  present,  then  2,  3 or  4 of  them  would  he  occupied:  we  need  to  know 
P(e|2£A£U),  where  E is  the  event  that  the  given  cell  is  occupied.  In  general, 
we  need  to  know  P ( E | a^^Aj^a^ ) , • We  have 

P(E|a^£A^a2)  = P(E  and  a^^A£a2)/P(a^£A£a2) 

ap  a 

= E P(E  and  A=a)  /E  P(A=a) 


= E'^  P(E|A=a)P{A=a)  / E^  P(A=a)  . 


For  this  we  need  the  marginal  distribution  of  A,  given  in  Table  3. 


[TABLE  3] 


2.3.  Example 


Let  us  now  apply  the  formula  to  the  configuration  of  Figure  1.  Receill  from 
the  Introduction  that  for  this  configuration  the  method  of  [3]  gave  - log^^ P = 11. 1+. 

To  make  the  adjustment  for  border  cells,  p^(a)  = P(Possibility  ilA=a) 
is  replaced  by  the  relevant  probability  of  the  form  P(Possibility  i|a,<A<a„). 

We  insert  the  relevant  cell  probabilities  on  Figure  2,  a copy  of  Figure  1. 

For  this  we  need  the  following  numbers,  obtained  from  the  tables: 
for  ending  ridges: — 

.356  P(e|A=0)  = .356(.199)  = -0708 

.356  P(e|0£A<,2)  = .356(.26i)  = .0929 


for  forks: — 


.163  p(e|a=o)  = .i63(.199)  = .032U 
.163  P(e|A=1)  = .163(.291)  = .0U7U 


for  empty  cells:  — 


1 - P(e|0£A<1|) 
1 - P(e10<A<3) 
1 - p(e|o<a<i) 

1 - P(e|1j<A<2) 
1 - P(e11<A<3) 


1 - .280 
1 - .2T3 

1 - .239 

1 - .309 
1 - .322 


.720 

.72'' 

.761 

.691 

.678 


[FIGURE  2] 


The  estimated  probability  is 

P = . 0321+^. 01+7':''’. 07 08^. 0929^. 667^. 678^. 691^. 7 09^^. 720^. 727 ^.761^. 801^  ; 

the  negative  log  of  this  is  - log^^ P = 12.5  . Compare  this  with  the  figure  of 
11.1+  given  by  the  approximation  based  on  an  assumption  of  independence  between 
cells.  The  difference  in  logarithms  is  1.1;  the  ratio  of  the  two  estimates  is 
thus  12.6.  This  difference  is  unimportant  since  we  are  interested  only  in  order 
of  magnitude.  Note  further  that  the  estimate  based  on  independence  is  a larger 
I probability,  i.e.,  it  is  conservative,  in  the  sense  of  giving  the  suspect  the 

benefit  of  the  doubt.  In  general,  independence  gives  too  much  weight  (too  low 
a probability)  to  configurations  with  a lot  of  clustering  of  occurrences.  In 
the  configuration  of  Figure  1 there  is  some  but  not  a great  deal  of  clustering. 

3.  DISCUSSION 


3.1.  Upper  Bounds 

Often  we  are  concerned  primarily  with  upper  bounds,  as  these  are  conservative^ 
i.e.,  they  err  in  favor  of  the  suspect.  An  upper  bound  on  the  probability  that  is 
highly  conservative  within  the  context  of  the  present  model  can  be  obtained  by 
multiplying  the  relative  frequency  of  each  occiirrence  by  .7ll*»  the  occupancy 
probability  for  1+  adjacencies,  and  replacing  the  probabilities  for  empty  cells  by 
.801,  the  probability  a cell  is  empty  given  that  all  its  neighbors  are  empty. 

This  gives  the  upper  bound 
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prints.  A print-dependent  notric  might  be  more  appropriate;  e.g.,  the  cell  size 
rouDd  be  tak€n  equal  to,  say,  the  average  distance  to  cross  six  ridge  lines. 

This  would  be  comparatively  di’-’icult  ;o  vcrk  with.  Use  of  the  model  of  the 
present  paper,  conditioning  on  occurrences  in  adjacent  cells,  provides  at  least 
a partial  control  for  variation  in  density  of  occurrences  associated  with 
variation  in  finger  sizes. 

APPIENDIX  A Data  Analysis  of  Dependence  between  Cells 

Square  blocks  of  9 cells,  3 cells  by  3 cells,  were  examined  to  determine 
the  extent  of  inter-cell  dependence.  The  data  set  of  [3]  yielded  8U5  separate 
such  blocks  of  cells.  For  i = 1,  2,  ...  , 8U5  blocks,  let  the  variable  y^  = 

1 or  0 according  as  the  center  cell  of  th  i-th  block  is  occupied  or  not,  and  let 
x^  be  the  number  of  adjacent  cells  which  are  occupied;  x^  is  between  0 and  8. 

Table  5 gives  the  cross-tabulatiOn  of  y and  x.  Table  6 gives,  for  each  value  of  x, 
the  proportion  of  y's  that  are  equal  to  1,  i.e.,  the  proportion  of  center  cells 
which  are  occupied. 

[INSERT  TABLES  5 AND  6 HERE.] 

The  probability  of  occupancy  increases  monotonically  with  x.  Such  absolute 
consistency  was  not  expected,  firstly  because  it  seons  so  rare  in  data  analyses 
and  secondly  because  it  was  thought  that  occurrences  in  most  of  the  adjacent  cells 
might  crowd  out  occurrence  in  the  center  cell. 

Table  7,  based  on  Table  5,  gives  observed  and  expected  values  of  the  number 
of  blocks  with  center  cells  occupied. 

[TABLE  7] 

The  value  of  the  chi-square  statistic  for  testing  independence  based  on 
Table  7 is  18.77  ( 6d.f.,  P<  .005  ).  The  decomposition  of  this  overall  value 

based  on  the  value  O.lU  of  the  correlation  coefficient  between  x and  y is  given 
in  Table  8. 


[TABLE  8] 
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APPENDIX  B Some  Remarks  on  Two-Dimensional  Processes 

The  idea  of  order  of  a two-dimensional  process  that  we  use  is  essentially  that 
discussed  by  P.K.  Bhattacharya  [2].  Consider  a stochastic  process  with  a two-dimensional 
parameter  (s,t).  In  order  to  calculate  probabilities  we  shall  need  to  define  a 
linear  ordering  on  {(s,t)}.  For  simplicity  assume  the  array  {(s,t)}  is  rectangular: 
s = 1,  2,  S,  t = 1,  2,  T.  We  shall  use  the  ordering  (l,l),  (l,2),  (l,T), 

(2.1) ,  (2,2),  ...,  (S,T),  which  might  be  called  "English-language  ordering." 

[Bhattacharya  uses  the  "lexicographic  ordering,"  (l,l),  (2,l),  (S,l),  (l,2) , 

(2.2) ,  (S,T).]  We  shall  say  that  a two-dimensional  process  is  of  order  r 

if  the  conditional  distribution  of  X , , ,,  given  all  the  random  variables  preceding 

s ,t 

it  in  the  ordering,  depends  only  upon  those  random  variables  X such  that 

S w 

2 2 2 

(s  - s')  + (t  - t')  £ r . (in  some  problems  it  may  prove  meaningful  to  use  a 

metric  different  from  Euclidean  distance.)  Thus  in  a first-order  process,  X depends 

sx> 

1/2 

upon  X and  X . In  a process  of  order  2 » X depends  upon  the  four 

S)  So 

variables  X ^,X  ,X  + . t > and  X . this  is  the  order  of  process 
employed  in  the  present  paper. 


APPENDIX  C Infinitely  Divisible  Random  Vectors 

Let  Y = (Y,  , Y , . . . , Y, „ ) be  the  observation  in  the  c-th  cell,  where, 

c Ic  2c’  10, c 

for  V = 1,  2,  . . . , 10  Galton  characteristics,  Y = number  of  occurrences  of  the 

VC 

v-th  characteristic  in  the  c-th  cell.  The  family  of  random  variables  {Y  , c=l,2,...,t 

c 

cells)  may  be  considered  as  a stochastic  process  in  the  plane,  a multivariate  point 
process.  It  is  reasonable  to  model  this  process  as  an  infinitely  divisible  process. 
Such  processes  have  been  studied  under  AFOSR  Grant  76-3050  (see  [^4])  and  are  being 
studied  under  AFOSR  Grant  77-3^5^-  Now  underway  is  a study  of  the  occurrences  of 
fingerprint  characteristics  as  a multivariate  Poisson  process,  a special  case  of 
an  infinitely  divisible  process.  Upon  completion  of  that  study,  infinitely  divisible 
processes  will  be  further  studied  and  will  be  applied  to  other  pattern  recognition 
problems  and  to  other  phenomena,  as  outlined  in  the  proposal  for  AFOSR  Grant  77-3^5**. 
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1.  Estimates  of  Prc^'Etiilities  and 
Relative  Frequencies  of  the  Possibilities 


Estimated  probability.  Estimated 
i p.  relative  frequency, 

^ Pi/d-Po^  = Pi/-23‘* 


0 

Einpty  cell 

.766 

1 

Island 

.0177 

.0756 

2 

Bridge 

.0122 

.0521 

3 

Spur 

.0071*5 

.0318 

1* 

Dot 

.0151 

.061*5 

5 

Ending  ridge 

.0832 

.356 

6 

Fork 

.0382 

.163 

T 

Lake 

.0061*0 

.0271* 

8 

Trifurcation 

.000582 

.0021*9 

9 

Double  bifur- 
cation 

.0011*0 

.00598 

10 

Delta 

.00198 

- 

.0081*6 

11 

Broken  ridge 

.0139 

.0591* 

12 

Other  multiple 
occurrence 

.0355 

.152 

1 

.0 

1.0 

1? 


2.  Number  of  Adjacencies  and  Probability  of  Occupancy 


\ 


Number  of  adjacencies 

Probability  of  occupancy 

0 

.199 

1 

.291 

2 

.333 

3 

.Uoo 

k 

.711* 

U.  Estimates  for  Border  Cells: 

\ 

Estimates  of  Probability  of  Occupancy, 
given  Partial  Information  about  Adjacent  Cells 


P(E|a^£A<a2) 

^2 

0 

1 

2 3 

k 

0 

.199  .239 

.261 

.273 

.280 

1 

.291 

.309 

.322 

.331 

2 

.333 

.352 

.369 

3 

.Uoo 

.I4UT 

!+ 

• Til* 

I 

1 


\ 
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I 

I 5.  Cross-Tabu  ation  of 

' Occupancy  of  Center  Cell 

I and  Number  of  Adjacent  Cells  Occupied 


X 

y 

0 

1 

Total 

0 

152 

28 

180 

1 

ITO 

1*5 

215 

2 

163 

1*5 

208 

3 

97 

29 

126 

U 

hh 

23 

67 

5 

23 

13 

36 

6 

7 

5 

12 

T 

0 

0 

0 

8 

0 

1 

1 

656 

189 

81*5  block 

y = 

1 if  given  (center) 

cell 

is  occupied 

' = 0 if  it  is  empty 

X = number  of  adjacencies  (number  of  adjacent  cells 
I occupied) 


~2‘.~ 


Probability  of  Occ  pancy  an  a Function 
of  Number  of  Adjnct  .t  Cells  Occupied 


Percent  of  blocks 
with  center  cells 

occupied:  15-6  20.9  21.6  23.0  3l*.3  36.1  hi. 7 100  22. h 


X = number  of  adjacencies 


T • Observed  and  Expected  Values 
of  Number  of  Blocks  with  Center  Cells  Occupied 


Observed  (Expected) 


X 

y 

0 

T 

X 

Total 

0 

152  (139.7) 

28  (I4O.3) 

180 

1 

170  (166.9) 

1^5  (1^8. 1) 

215 

2 

163  (161.5) 

1*5  (I46.5) 

208 

3 

97  (97.8) 

29  (28.2) 

126 

U 

(52.0) 

23  (15.0) 

67 

5 

23  (27.9) 

13  ( 8.1) 

36 

6 or  more 

7 (10.1) 

6 ( 2.9) 

13 

8.  Decomposition  of  Chi-Square 
according  to  Correlation  between  x and  y 


Source  of 

variation 

d.f . 

Value  of 

chi-square 

Overall  6 

18.77 

(P<.005) 

nr^  1 

16.65 

(P<.005) 

Residual  5 

2.12 

( .80<P<.85) 

Residual 


1 2 3 4 5 6 

A 0 0 0 0 0 0 
B E 0 0 O E 0 
C 0 0 F 0 0 0 
D 0 0 0 0 0 0 
E O 0 0 E 0 0 
F -0  0 0 0 F 0 
C 0 E 0 0 O 0 
H 0 


Figure  1«  Conf iguretlon  of  43  cells  with  4 ending  ridges 
end  2 forks.  0 ■ empty  cell,  E *>  ending  ridge, 
F ■ fork. 


-2('- 


1 

2 

3 

U 

5 

6 

A 

.720 

.727 

.727 

.727' 

.727 

.727 

B 

.0929 

.709 

.801 

.801 

.0708 

.691 

C 

.678 

.709 

.0324 

. 667 

.709 

.691 

D 

.691 

.709 

.709 

.709 

.801 

.761 

E 

.691 

.801 

.801 

.0708 

.709 

.761 

F 

.691 

.801 

.709 

.709 

.0474 

.691 

G 

.691 

.0708 

.709 

.709 

.709 

.691 

H .691 


Figure  2.  Cell  probabilities  for  configuration 
of  Figure  1, 
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